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Constraining nucleon high momentum in nuclei 
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Recent studies at Jefferson Lab show that there are a certain proportion of nucleons in nuclei 
have momenta greater than the so-called nuclear Fermi momentum pf. Based on the transport 
model of nucleus-nucleus collisions at intermediate energies, nucleon high momentum caused by the 
neutron-proton short-range correlations in nuclei is constrained by comparing with n and photon 
experimental data and considering some uncertainties. The high momentum cutoff value p m ax < 
2 pf is obtained. 

PACS numbers: 25.70.-z 


The picture of nucleons have maximal momentum - 
so called the Fermi momentum pp — in a nuclear system 
and roughly move independently in the mean field cre¬ 
ated by their mutually attractive interactions has been 
established since the 1950s. However, recent proton- 
removal experiments using electron beams with energies 
of several hundred MeV showed that only about 80% nu- 
cleonsparticipate in this type of independent particle mo¬ 
tion |lH3j. And high-momentum transfer measurements 
have shown that nucleons in nuclei can form pairs with 
larger relative momenta and smaller center-of-mass mo¬ 
menta 0,0- This is interpreted by the nucleon-nucleon 
tensor interaction in short range 00- The nucleon- 
nucleon short-range correlations (SRC) in nuclei leads to 
a high-momentum tail in single-nucleon momentum dis¬ 
tribution above 300 MeV/c [5412} . And interestingly, the 
high-momentum tail’s shape caused by two-nucleon SRC 
is almost identical for all nuclei from deuteron to very 
heavier nuclei [TaI - ITgI . i.e., roughly exhibits a C/k 4 tail 
fl7l - [20l |. Nucleon momentum distributions at even higher 
momenta are due to three or many-nucleon correlations. 
This part of momentum-distribution probability was de¬ 
duced to be less than 1% (2lj . We thus in this study 
neglect this kind of high-momentum nucleons caused by 
many-nucleon short-range correlations. 

In the high-momentum tail (HMT) of nucleon momen¬ 
tum distribution, nucleon component is strongly isospin- 
dependent, i.e., the number of n-p SRC pairs is about 
18 times that of the p-p or n-n SRC pairs 13), thus in 
neutron-rich heavy nuclei proton has greater probability 
than neutron to have momenta greater than the nuclear 
Fermi momentum |20j. In neutron stars, the number of 
protons only has a small proportion. The above n-p SRC 
in neutron stars will cause proton average kinetic energy 
far greater than neutron’s [22| . And the stronger the n-p 
SRC is, the larger the difference of proton and neutron 
average kinetic energy is seen. 

Nucleon spectral function provides fundamental infor¬ 
mation on the dynamics of nucleon in nuclear medium. 
The nuclear momentum distribution can be obtained 
from the spectral function by integrating over the excita¬ 
tion energy fl3l [23}. The high-momentum nucleons come 
predominantly from the high excitation energy regime 
of the spectral function. The high-momentum cutoff pa¬ 
rameter A (= Pmax/PFi he., the ratio of nucleon maximal 


momentum over the nuclear Fermi momentum) was first 
introduced by Hen et al. as a free parameter in the Cor¬ 
related Fermi Gas model an analytical approximation 
for the momentum distribution of nucleon in symmetric 
nuclei and nuclear matter jl8} , to avoid divergence when 
calculating nucleon average kinetic energy assuming a 
C/k 4 dependence for the high-momentum tail. There¬ 
fore, the implication of the value of this effective param¬ 
eter A is the determination of the average kinetic energy 
of nucleons. 

The value of average kinetic energy of nucleons and 
the high-momentum tail of nucleon distribution in nuclei 
surely affect the yields of 7r, K , p and nucleon emission in 
heavy-ion collisions at intermediate energies. The isospin 
dependence of nucleon high-momentum distribution def¬ 
initely affects transport calculations of the symmetry- 
energy sensitive observables. More specifically, a low 
(high) value of average kinetic energy of nucleons causes 
a small (large) number of meson production in trans¬ 
port calculations owing to low (high) value of collision 
energy of nucleon pairs. Because in neutron-rich heavy 
nuclei protons have greater probability than neutrons to 
have momenta greater than the nuclear Fermi momen¬ 
tum, the high-momentum tail of nucleon momentum dis¬ 
tribution affects values of n~/p + ratio and the differ¬ 
ence of neutron and proton elliptic flows [24} . Values of 
average kinetic energy of neutrons and protons in nu¬ 
clei also strongly affect nuclear kinetic symmetry energy 
EQ1, the latter is known pla ys crucial role in both 
nuclear physics and astrophysics [25} . 

In fact, one can deduce the high-momentum cutoff 
parameter A from nucleon momentum distribution in 
deuteron [l7|,[13 or from the high-energy electron scatter¬ 
ing measurements [2C| [2l| . But the high-energy electron 
scattering measurements main ly p robe the nucleon mo¬ 
menta at the surface of nuclei [26} . Since the production 
of 7r + meson in nucleus-nucleus collisions at intermedi¬ 
ate energies is mainly from proton-proton collision [27} . 
and energetic neutron-proton scattering produces hard 
bremsstrahlung photon [28l - [30} . the high momentum of 
nucleon in projectile or target should affect 7r + and hard 
photon productions. In this study, we use hadronic 
probe 7r + meson and electromagnetic probe hard pho¬ 
ton in nucleus-nucleus collisions to probe nucleon high- 
momentum cutoff value. 
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To obtain the high-momentum cutoff value of nu¬ 
cleon in nuclei by nucleus-nucleus collisions, we use the 
Boltzmann-Uehling-Uhlenbeck (BUU) transport model 
a which has been very successful in studying heavy- 
ion collisions at intermediate energies. The BUU trans¬ 
port model describes time evolution of the single particle 
phase space distribution function f(f,p, f), which reads 

df 

+ VpE ■ - V P E ■ Vpf = I c . ( 1 ) 

The phase space distribution function f(r,p,t ) denotes 
the probability of finding a particle at time t with mo¬ 
mentum p at position r. The left-hand side of Eq. m 
denotes the time evolution of the particle phase space 
distribution function due to its transport and mean field, 
and the right-hand side collision item I c accounts for the 
modification of phase space distribution function by elas¬ 
tic and inelastic two body collisions j3ll - l33j . E is a par¬ 
ticle’s total energy, which is equal to kinetic energy Eki n 
plus its average potential energy U. While the mean-field 
potential U of single particle depends on its position and 
momentum of the particle and is given self-consistently 
by its phase space distribution function f(r,p,t). 
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FIG. 1: (Color online) Momentum distribution n(k) of neu¬ 
tron and proton in nucleus yg 7 Au with normalization condi¬ 
tion 4nf 0 Xk r n p ,n(k)k 2 dk = 1 as well as 47r J/.p F n p (k)k 2 dk 

~ 25% and 4 - 7 T f X n F n n (k)k 2 dk ~ 17%. 


In the used BUU model, nucleon spatial distribution 
in initial colliding nuclei is given by [311 ] 

r = R(x \) 1 ^ 3 ; cosQ = 1 — 2x2', <t> — 2nx3. (2) 

x = rsinOcoscj)', y = rsin9sin(f>; z = rcosO. (3) 

Where R is the radius of nuclei, X\,X 2 ,x 3 are three in¬ 
dependent random numbers. Since there is a rough 20% 


depletion of nucleon momentum distribution inside the 
Fermi sea 3, EH H3, [H a • we let nucleon momen¬ 
tum distributions in the high-momentum tail n HMT (k) ~ 
1 /k 4 [T?! and J^ F n HMT (k)k 2 dk/ f^ kp n(k)k 2 dk ~ 20 % 
and keeping rip MT (k)/n^ MT (k ) ~ N/Z (N and Z being 
the neutron and proton numbers of a nucleus ) @, Eiil EH ■ 
As shown in Fig. [0 there are about 25% (17%) protons 
(neutrons) with momenta larger than the proton (neu¬ 
tron) Fermi momentum. With this nucleon momentum 
distribution, the average kinetic energy of nucleons in this 
study increases roughly several MeV comparing to that 
with ideal Fermi-Gas model. We thus neglect this dif¬ 
ference in heavy-ion collisions at 400 MeV/nucleon beam 
energy. 

In this model, an isospin- and momentum-dependent 
mean-field single nucleon potential is used j3f| 123 ] , which 
reads 


U(p,8,p,r) = 


A u {x)— + Ai{x) — 

Po Po 

— 1 — xS 2 ) — 8xt —-——- 5p T ' 

Po o + 1 Po 


+ 


2CV,t /' ,3 j fr(r,p) 


PO 


d 3 p - 


1 + ip-p') 2 / A 2 


+ f d 3 p' /r ' ( ^2..„ , (4) 


Po 


1 + (p-p') 2 /A 2 


where t,t' = l/ 2 (— 1 / 2 ) for neutrons (protons), S = 
( p n — Pp)/{p n + Pp) is the isospin asymmetry, and p n , p p 
denote neutron and proton densities, respectively. The 
parameter values A u (x) = 33.037 - 125.34a; MeV, Ai(x) 
= -166.963 -(- 125.34a; MeV, B = 141.96 MeV, C TyT = 
18.177 MeV, C T>T > = -178.365 MeV, a = 1.265,’and 
A = 630.24 MeV/c are obtained by fitting empirical con¬ 
straints of the saturation density po = 0.16 fm~ 3 , the 
binding energy Eo = -16 MeV, the incompressibility I\q 
= 230 MeV, the isoscalar effective mass m* = 0.7 to, the 
single-particle potential = 75 MeV at infinitely large 
nucleon momentum at saturation density in symmetric 
nuclear matter, the symmetry energy S(p) = 30 MeV (we 
let kinetic symmetry energy roughly be 0 MeV 38]) and 
the symmetry potential U^ m = -100 MeV at infinitely 
large nucleon momentum at saturation density. f T {f,p) 
is the phase-space distribution function at coordinate r 
and momentum p and solved by using the test-particle 
method numerically. Different symmetry energy’s stiff¬ 
ness parameters x can be used in the above single nu¬ 
cleon potential to mimic different forms of the symme¬ 
try energy predicted by various many-body theories j39] 
without changing any property of the symmetric nuclear 
matter and the symmetry energy at normal density. In 
this study, however, both 7 r + production in Au + Au col¬ 
lisions and hard photon production in C + C collisions 
are in fact not sensitive to the symmetry energy param¬ 
eter x. 

According to baryon effective mass, the isospin- 
dependent baryon-baryon ( BB ) scattering cross section 
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in medium j s reduced compared with their free- 

space value a BB e by a factor of [33 ] 


Rmedium(p,5,p) = / (T*™ 

= (Pbb/pbb) 2 , (5) 


where /is b and p* BB are the reduced masses of the collid¬ 
ing baryon-pair in free space and medium, respectively. 
This form of reduced elastic baryon-baryon scattering 
cross section in medium agrees well with our recent study 
[2E|. Since the inelastic baryon-baryon scattering cross 
section in medium is less known but crucial for 7r produc¬ 
tion [io} . we in the present model extend the above re¬ 
duced factor RmeAinm.(Pi d, p) to inelastic baryon-baryon 
scattering cross section [4l| . Other treatments related to 
7 T production are similar to that in Ref. [42}. 

For hard photon production from neutron-proton 
bremsstrahlung, we use the prediction of the one boson 
exchange model by Gan et al. [28l - f3(i [43] 


p 7 = 


dN _ 7 _ 6 (1 - y 2 )“ 


ds~, 


= 2.1 x 10” 


( 6 ) 


where y = £j/E max , a = 0.7319 — 0.5898/%, e 7 is energy 
of emitting photon, E max is the energy available in the 
center of mass of the colliding proton-neutron pairs, /% 
is the initial velocity of the proton in the proton-neutron 
center of mass frame. The pauli-blockings of final state 
scattering neutron and proton in pn —» pn'y process are 
taken into account [44| . 
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FIG. 2: (Color online) The number of produced 7r + meson as 
a function of high-momentum cutoff parameter A in the Au 
+ Au collisions at, respectively, 0.4 and 1 GeV/nucleon beam 
energies. 

Fig. [2] shows 7T + production as a function of high- 
momentum cutoff parameter A of colliding nuclei in the 
Au + Au collisions at 0.4 and 1 GeV/nucleon incident 
beam energies, respectively. One can clearly see that as 
the high-momentum cutoff parameter A increases, more 
7 r + ’s are produced. Larger high-momentum cutoff pa¬ 
rameter A causes larger nucleon average kinetic energy, 
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especially proton average kinetic energy [2Cj, thus the 
average center-of-mass energy of proton-proton collision 
also becomes larger. As a consequence more 7r + ’s are 
produced in nucleus-nucleus collision [27j]. This is the 
reason why one sees in the upper panel of Fig. [2] more 
7 r + ’s are produced with large high-momentum cutoff pa¬ 
rameter A. As incident beam energy increases, the initial 
movement of nucleons in nuclei becomes less important 
in nucleus-nucleus collisions. We thus see, in the lower 
panel of Fig. [2l at 1 GeV/nucleon incident beam energy, 
7 r + production is less sensitive to the high-momentum 
cutoff parameter A (At 0.4 GeV/nucleon, the sensitivity 
of 7r + production to A is about 10 times that of 7r + at 
1 GeV/nucleon). Fig. [2] shows A < 2 is favored by the 
FOPI data 0. 



FIG. 3: (Color online) Constraints on the high-momentum 
cutoff parameter A by the ratio of 7r + productions in the Au 
+ Au collision at low and high beam energies. The shadow 
region denotes the FOPI data [45} . 


Since 7r + production is sensitive to the high- 
momentum cutoff parameter A only at relative low beam 
energy, one can construct the ratio of 7r + ’s produced from 
low and high incident beam energies to probe the high- 
momentum cutoff parameter. This ratio is expected to 
reduce the system errors in some degree. Shown in Fig. [3] 
is the ratio of 7r + productions at 0.4 and 1 GeV/nucleon 
incident beam energies as a function of high-momentum 
cutoff parameter A. As expected, the ratio of 7r + multi¬ 
plicities produced respectively from low and high beam 
energies are still very sensitive to the high-momentum 
cutoff parameter A. By comparison with the FOPI pion 
production data [45}, a constraint of 1.5 < A < 2.5 is 
obtained. This result is surprisingly similar to that in 
Ref. [H|. 


The hadronic probe 7r + production inevitably suffers 
from distortions due to the strong interactions in the fi¬ 
nal state. Ideally one would like to have more clean ways 
to probe the high-momentum tail of nucleon in nucleus. 
Photons interact with nucleons only electromagnetically, 
once produced they escape almost freely from the nu¬ 
clear environment in nuclear reactions. In this regard, we 
also use hard photon production to constrain the high- 
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momentum cutoff parameter A. 

Hard photon production in heavy-ion reactions at 
beam energies below 200 MeV/nucleon had been in fact 
extensively studied both experimentally and theoretically 
[28l - l3ll . I46l - l48| . The TAPS collaboration carried out a 
series of comprehensive measurements studying in detail 
the properties of hard photons [49l - [52| . Theoretically, it 
was concluded that the neutron-proton bremsstrahlungs 
in the early stage of the reaction are the main source of 
high energy 7 rays jH, [54|. And it was demonstrated 
that the hard photons can be used to probe the reac¬ 
tion dynamics leading to the formation of dense matter 
[44 l55l - [58j . And effects of the nuclear Equation of State 
(EOS) on the hard photon production were found small 
[59j . Although the input element ary pn —> pn"f probabil¬ 
ity is still model dependent [43l . [60l463j , the experimen¬ 
tal data can be described reasonably well theoretically 
within a factor of 2 [48}. And the experimental efforts 
have the potential to improve the situation significantly 
in the near future [64}. 



FIG. 4: (Color online) Inclusive photon production cross 
sections (e 7 > 150 MeV) in 12 C+ 12 C collisions at the beam 
energy of 60 MeV/nucleon. The symbols stand for BUU cal¬ 
culations with, respectively, A = 1 (i.e., without HMT), 1.5, 
2, 2.5. The shadow region denotes experimental data [2a. [4611 . 

Fig. [4] shows comparison of theoretical inclusive hard 
photon production cross sections in 12 C+ 12 C collisions 
and the experimental data [28|, [46|. Since the hadronic 
probe 7 r + production constrained the high-momentum 
cutoff parameter A between 1.5 ~ 2.5, we use A = 1.5, 
2, 2.5 as the BUU calculations with the high-momentum 
tail. As comparison, we also calculated the case with¬ 
out high-momentum tail (A = 1). From Fig. 21 it is seen 
that the hard photon production cross section in heavv- 
ion collisions is very sensitive to the HMT of nuclei [30| ■ 
The BUU calculations with A = 2 and 2.5 are larger than 
the experimental hard photon production cross section. 


And the case of BUU calculation without HMT is some¬ 
what lower than the experimental data. From Fig. [H 
it is seen that the electromagnetic probe hard photon 
constrains the the high-momentum cutoff parameter to 
be A < 2. Combining the constraints from hadronic 
probe 7 r + production (shown in Fig. [3]) and that from 
electromagnetic probe hard photon production (shown 
in Fig. SI), we can conservatively conclude that the value 
of the high-momentum cutoff parameter A in nuclei is less 
than 2.5 and the overlap-area is A < 2, which is smaller 
than that deduced in other lecture [l7i - [l9l| . 

A small value of lambda implies lower average nucleon 
kinetic energy. The lower average nucleon kinetic energy 
implies smaller collision energy of nucleon pairs in trans¬ 
port calculations. This causes small number of meson 
production in heavy-ion collisions at low and intermedi¬ 
ate energies. And it also cause smaller number of ener¬ 
getic nucleon or meson emissions, a high value of n~/ ti + 
ratio [24l ] and small number of hard photon production 
in heavy-ion collisions at low and intermediate energies. 
A small value of lambda also implies a relatively larger 
nuclear kinetic symmetry energy, thus causes the reduc¬ 
tion of nuclear symmetry potential [ 18[ |. The reduction 
of nuclear symmetry potential in heavy-ion collisions at 
intermediate energies decreases the sensitivity of isospin- 
sensitive observables. 

Because the high-momentum tail of nucleon momen¬ 
tum distribution is in fact caused by the short-range cor¬ 
relations of nucleons, while in our transport model (be¬ 
sides nucleon-nucleon or nucleon-meson collisions and nu¬ 
clear pauli-blockings) only a mean-field potential is used. 
Thus the nuclei in the evolution before collision may be 
instable. Lacking of the binding of high-momentum nu¬ 
cleons in nuclei, the shape of initial distribution of nucle¬ 
ons in momentum space may be changed and energetic 
nucleons may escape out of the nuclei [4l| ■ Therefore, 
the average kinetic energy of nucleons in the reaction 
system decreases and then cause less 7 r + meson or hard 
photon productions. The increased stability of collid¬ 
ing nuclei may cause somewhat more 7 r + meson or hard 
photon productions (less emission of high-momentum nu¬ 
cleons corresponds to a higher average kinetic energy of 
nucleons in nuclei). Thus considering stability of collid¬ 
ing nuclei, more pion and hard photon may be produced 
than that in the present work. In a word, after consid¬ 
ering stability of colliding nuclei and comparing to the 
same experimental data, the high-momentum cutoff pa¬ 
rameter A should be smaller than our current constraint. 
Also uncertainty of the mechanism of hard photon pro¬ 
ductions may affect the conclusion here. The probability 
of hard photon production from the semiclassical hard 
sphere collision model [3ll l47i . l48j | will give somewhat 
more photon production [28[| , thus also require a smaller 
A value to explain the experimental data. Furthermore, 
the off-shell transport of particle production in the BUU 
model may also cause more 7 r + meson and hard photon 
productions [6514671 , thus also require a smaller A value to 
fit the experimental data. As in this work the A value is 
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constrained to be A < 2, while considering all the above 
factors, the A parameter should be not larger than 2. 

In conclusions, based on the nuclear transport model, 
we studied how the high-momentum cutoff parameter A 
affects 7T and hard photon productions in nucleus-nucleus 
collisions at intermediate energies. It is found that 7r + 
and hard photon productions in nucleus-nucleus collision 
at lower beam energy is very sensitive to the value of the 
high-momentum cutoff parameter A. By comparing the 
BUU’s 7r + and hard photon productions with experimen¬ 
tal data and considering some uncertainties, a constraint 
of high-momentum cutoff value A < 2 is obtained. 

Constraints on the high-momentum cutoff parameter A 


in nuclei have implications in the studies of nuclear force 
at short distance, in the construction of nuclear transport 
model of heavy-ion collisions at intermediate energies, in 
the studies of equation of state of dense nuclear mat¬ 
ter and the nuclear symmetry energy at suprasaturation 
densities or in the study of the physics in neutron stars, 
etc. 

The work was carried out at National Supercomputer 
Center in Tianjin, and the calculations were performed 
on TianHe-lA. The work is supported by the National 
Natural Science Foundation of China under Grant Nos. 
11375239, 11435014. 
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